Fokker - Planck equation in curvilinear coordinates. Part 2 



Igor A. Tanski 
Moscow, Russia 
povorot2 @ infoline. su 



V&T 



ABSTRACT 



The aim of this paper is to derive Fokker - Planck equation in curvilinear coordi- 
nates using physical argumentation. We get the same result, as in our previous article [1], 
but for broader class of arbitrary holonomic mechanical systems. 



1. Lagrange equation 

In our previous article [1] we derived Fokker - Planck equation in curvilinear coordinates from its 
Cartesian form using transformation of space coordinates and velocities. This way does not show physical 
sense of final equation and involved variables. 

We begin from Lagrange equation [2] for holonomic mechanical system 



where T - kinetic energy of the system; 
q' - generalized coordinates; 
Qi - generalized force. 

In the following, we use the Einstein summation notation for implicit sums. 

We have following expression for kinetic energy 



Inertia tensor gy plays here the same part as the metric tensor in our article [1]. It depends on gener- 
alized coordinates and transforms as covariant tensor by transformations of generalized coordinates q l . We 
use this tensor to lower indices, that is to build one to one correspondence between contravariant and 
covariant tensors. If the considered mechanical system is a particle of unit mass, both tensors coincide. We 
therefore use in present article the same notation gy for inertia tensor as in article [1] for metric tensor. 
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Velocities q' transform as components of contravariant vector. 
Momentum vector 

dT 

Pl = w = 8u q J . (3) 

transforms as covariant vector and generalized force Q t also. 

The determinant g = det(gy) defines 2 measures of volume: measure of volume in configurations 
space ~{gdq l dq 2 ■ ■ ■ and measure of volume in velocities space Jgdv l dv 2 ■ ■ ■. Both measures are invariants 
of arbitrary transformations of generalized coordinates. 

The product of these measures is measure of volume in phase space. According to Liouville theorem 
it is invariant of infinitesimal transformation, induced by motion of conservative system. 

Let us represent generalized force as a sum of random force and friction force, which is propor- 
tional to velocity 

0, = R f - a 8ij qj. (4) 

Lagrange equation (1) takes the form 

d ( j)_ 1 dgpq _ 



*^J-2-^ W = *'- a *^- (5) 

This differential equation of the second order we can write down as a system of two differential equa- 
tions of the first order. We can do it in many ways, depending on our choice of unknown variable. 2 ways 
are considered below: contravariant and covariant components of velocities. 



2. Contravariant components of velocities 

Equation (5) with contravariant components of velocities as unknown variables is: 

Expand the first term in (6) 

d ( A d 8y i , dyJ d Sij p i , dv J 

jt ^ v + 8,j = v vv +8ij ^' (7) 

dV 

solve (7) for — and get 
at 
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dv J 
dt 



dqP 



(8) 



where g' J is inverse matrix 



S p 



(9) 



and S p is Cronecker delta . It is known [3], that g ij builds contravariant tensor. 



Terms with derivatives of inertia tensor are equal to Christoffel symbol [3] 



1 dgpq p a ^8 is n s 1 ^8pq n q 1 ^8qi p q 1 ^8ip n q r n q 

— ^ v' v q - - — v' v = — z-^- v' v q ■=- JL - v' v q ^r— - v' v q = -Vi, na V F V q . 

2 dq 1 dqP 2 dq' 2 dqP 2 dq q pq 



(10) 



In this way we obtain the full system with contravariant components of velocities as unknown vari- 



ables: 



dq' 
dt 



: V, 



(11) 



dv 



ih g ll R i -av 1 -V pq v p v'>. 



(12) 



If generalized force is zero, equations (11-12) state, that absolute derivative [3] of v 1 vector with 
Sv j 

respect to metric g u is zero = 0. This implies, that the paths of free movement of mechanical system 

8t 

are geodetics of metric g tj , defined by its inertia tensor. 

Let us consider now not a single system, but ensemble - a large set of noninteracting identical sys- 
tems. Let us denote h - density of systems in joint space of configuration coordinates and velocities . This 
means, that n" dq 1 dq 2 ■ ■ ■ dv 1 dv 2 is total number of systems in elementary cell, h satisfies evolutionary 
equation. 



dh d{hV) d 
dt dq' 3w 



= 0. 



(13) 



or 



^ + v> ^ + ^- {ng ij Ri V |^ [av> + T j pa v"v q )- 3an - 2hT\ v v p = 0. (14) 
dt dq 1 dviy ) dvJ { Pq J JP 

We calculate random force from the law 

dn 

nR^-k — . (15) 
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analogous to diffusion (Fick) and heat transfer (Fourier) laws. 

This phenomenological law we consider as definition of perfectly random force. 

We see, that in the case of zero forces (a = and k = 0) substitution n = g satisfies equation (16). 
This follows from the known identity [3] 

r> = - *H (17) 

This means, that flow without external forces preserves phase space volume. This is Liouville theo- 
rem for our case (we do not consider potential forces). 

These considerations suggest another definition of density. Let us define 

h = n det(gy) = ng. (18) 

We get for new density n equation 



dn j. dn , „ „ dn , dn ,, d 2 n 

di +v w- T > Pv9 M- av w- 3an = kg dVd^- (19) 



which coincides with equation (19) from our article [1]. 

Scaling (18) means, that new density n relates not to product of differentials, but to invariant product 
of volumes. 



3. Covariant components of velocities 

Covariant components of velocities are also components of momentum vector (see (3)). Therefore 
we use notation p l for them. 

Equation (5) takes the form 

dpi 1 dgpq .„.„ /orv\ 

This leads to following system of differential equations: 

it = 1 2 d §qfs pm 8 qn p m Pn + R l -a Pl , (21) 
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dq l 
dt 



: 8 Pn 



(22) 



We use the ensemble concept again. Let us denote n - density of systems in joint space of configura- 
tion coordinates and momentum vectors . This means, that n dq l dq 2 ■ ■ ■ dp x dp 2 • • • is total number of sys- 
tems in elementary cell. 

The measure of volume in velocities space is — dp\dp 2 therefore the measure of volume in 
phase space is dq 1 dq 2 ■ ■ ■ dp x dp 2 and n so defined coincides with density from the section 2. 



n satisfies following evolutionary equation: 



dn | d(ng im p m ) | d 
dt dq' dp i 



n\\ d ^ L g pm g qn p mPn + R i -ap i 



0. 



(23) 



Second term in (23) expands as follows: 



d(ng im Pm ) im dn dg'» 
— ^i— = 8 Pm^ + nP^ 



(24) 



Third term in (23) simplifies to 



j-fl dgpq pm qn "| = 1 ^gpq pj qn , l ^Spq pm qi _ ^Spq pj qn 

d Pi {2 dq t 8 8 PmPn ) 2 dqj 8 8 P " 2 d qi 8 8 Pm dqj 8 8 Pn ~ 



(25) 



So we have 



¥ + * W + ^ W + dp-, [2 if 8 8 PmP " ~ aP ' l + 



(26) 



+n ^ 8 Pi 8 qn Pn ~3an + ^- (n R f ) = 0. 
dqj dp. 



The sum of coefficients by n is equal to zero 



^1 + ^LeI g pj g qn _ o. 
dq' dqj 



(27) 



To prove this, we begin from identity 8'"8nk = <>k- Differentiate this by q' and obtain 



(28) 
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Then lift the index and get the desired result 



^ + /V"^f=0. (29) 



For the random force we have expression (15). It follows 



dn dn dp-, dn 

nR ^- k ^ = - k W j ^ = - k8 ^- (30) 



So we have 



dn . Am„ dn . dn ( l d 8 fl pm qn __ } .._ dn „ , d z n 



dt + ^ w + w t (2 it 8PmgqnpmPn r aPi *ir 3an= k8ij mr/ (31) 

Now it remains only one step to the final result: 

dn , dn „ „, dn dn , d 2 n 

di +s p -w k+ni8 PpP «wr aPk wr 3an = k8ik dpW k - (32) 

To prove equivalence of (31) and (32), we only need to check identity 

PmPn{\ ^ 8 Pm 8 q " ~ = 0. (33) 



2 dq 

It follows directly from definition of Christoffel symbol 



I ' = 1 fdgfa- | %g sk _ dg^ \ km , , 

2\dq s dq t dq k f 



This gives 



1 fdg k i dg is \ km 
- PmPn 2{d qs ~Wj 8 ( } 



But expression inside of parenthesis is antisymmetric on indices k, s and expression outside of paren- 
thesis is symmetric. Therefore the sum is zero, q.e.d. 
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DISCUSSION 

We proved, that equations (19) and (32) are forms of Fokker - Planck equation for contravariant and covari- 
ant components of velocities respectively. They coincide with forms, derived in our article [1]. 

We proved (19) and (32) not only for particle movement in curvilinear coordinates, but for broader class of 
arbitrary holonomic mechanical systems. 

Usually authors [4, 5] begin their studies from Hamilton's equations. This is not our way, because we wish 
to consider friction force and random force from the very beginning. Therefore we need instrument like 
Lagrange equations. 

Our equations (19) and (32) are analogous to Liouville equation [5] for density in phase space. This builds 
connection of Fokker - Planck equation approach with another works in the field of statistical mechanics . 
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